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Using Monte Carlo (MC) computer simulations we explore the self-assembly and ordering behavior 
of a hybrid, soft magnetic system consisting of small magnetic nano-spheres in a liquid-crystalline 
(LC) matrix. Inspired by recent experiments with colloidal rod-matrices we focus on conditions 
where the sphere and rod diameters are comparable. Already in the absence of a magnetic field, the 
nematic ordering of the LC can stabilize formation of magnetic chains along the nematic or smectic 
director, yielding a state with local (yet no macroscopic) magnetic order. The chains, in turn, 
increase the overall nematic order, reflecting the complex interplay of the structure formation of 
the two components. Increasing the sphere diameter the spontaneous uniaxial ordering is replaced 
by biaxial lamellar morphologies characterized by alternating layers of rods and magnetic chains 
oriented perpendicular to the rod’s director. These novel ordering scenarios at zero field suggest 
complex response of the resulting hybrid to external stimuli such as magnetic fields and shear forces. 


Mixtures of magnetic nanoparticles (MNP) in liquid- 
crystalline (LC) matrices attract attention in fundamen¬ 
tal physics and material science since about four decades 
ago: in 1970, Brochard and de Gennes [1] predicted, 
based on free energy considerations, various new effects 
including spontaneous magnetization in the liquid state, 
’’compensated” phases and giant field-induced effects. 
Whereas the field sensitivity has been investigated and 
verified in many experiments [2-8], new interest was stim¬ 
ulated by the recent discovery of spontaneous magnetic 
ordering in a hybrid system of large, micron-sized mag¬ 
netic plates embedded in a lyotropic nematic LC [9]. A 
crucial ingredient for the observed behavior are the LC- 
mediated interactions stemming from local distortions of 
the LC director via the presence of the plates, whose size 
is much larger than that of the LC molecules. Indeed, 
LC-mediated elastic interactions and their consequences 
for colloidal self-assembly and ordering have intensively 
studied by experiments [10] and simulations [11, 12] for 
many (non-magnetic) colloid/LC mixtures. 

In the present letter we focus on the opposite case of 
small MNPs in LC matrices, where the sizes are compa¬ 
rable and distortions of the LC do not play a significant 
role. An example of such systems are lyotropic suspen¬ 
sions of colloid pigment rods and magnetite MNPs which 
have been recently studied to investigate field-induced or¬ 
dering effects [13, 14]. However, although there is strong 
interest in using such hybrids as new, stimuli-responsive 
or even adaptive materials [9], the microscopic structure 
of small MNPs in LC matrices is essentially unexplored; 
although this is clearly a prerequisite in terms of advanc¬ 
ing new functionalities. Experimentally, the challenge is 
to actually ’’see” the relatively small MNPs with typ¬ 
ical diameters of 10-15 nm; a solution of this problem 
might be provided by recently developed tomography 
techniques [15]. Suprisingly, however, there is also a lack 
of computer simulation studies (contrary to the manifold 
of studies on non-magnetic rod-sphere mixtures [16-18]); 


indeed, to our knowledge ferronematics have so far only 
been studied on the basis of Landa-de Gennes free energy 
approaches [5, 19, 20]. Thus, microscopically, even basic 
questions are currently open, such as: do we find any 
self-assembly or collective ordering of the MNPs, how do 
they influence the LC and vice versa, and what is the 
role of the rod-sphere size ratio? 

To explore these questions we here employ canoni¬ 
cal MC simulations [21] of a simple binary model sys¬ 
tem consisting of N r rods, represented as prolate, non¬ 
magnetic ellipsoids and N s MNPs, represented as dipo¬ 
lar soft spheres (DSS), i.e., spheres with an embedded 
point dipole moment fi i ( i = 1,..., 7V S ) in their center. 
The rods interact via the anisotropic Gay-Berne (GB) 
potential involving both, repulsive and attractive (van- 
der-Waals) interactions, and we use a standard parame¬ 
terization [16]. The length (Z)-to- width (a) ratio is set 
to l/d = 3; this value is well-studied for pure GB sys¬ 
tems [22], and it is also in the experimentally accessible 
range for LC/MNP mixtures [14]. The dipolar spheres 
have diameter a s . We focus on systems with comparable 
sizes with cr* = cr s /cr = 1 and cr* = 2. The spheres inter¬ 
act via a combination of a soft-sphere potential and the 
conventional dipole-dipole potential. The reduced dipole 
moment is set to /i* = /a/ yeo of = 3 and we consider re¬ 
duced temperatures T* = k^T/e o (where k b is the Boltz¬ 
mann constant and eo is the soft-sphere energy parame¬ 
ter) in the range 0.6 to 2.5. Thus, the coupling parameter 
A = pi 2 /ks Ter g takes values between ^3.5 and ~ 15, con¬ 
sistent with values characterizing real, strongly coupled 
MNPs [23-25]. Long-range interactions are treated with 
the three-dimensional Ewald sum [26]. Finally, the in¬ 
teractions between rods and spheres are described via a 
modified GB potential [16, 27], parametrized such that 
the spheres favor the sides of the rods. We have exam¬ 
ined systems for a variety of total densities p* = 7Vcr 3 /U, 
compositions x a = N a /N (where a=r, s), and resulting 
partial volume fractions (f) a = N a v a /V (where v a is par- 
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tides’ and V the total volume). For all chosen systems, 
we have performed cooling and heating series by varying 
the temperature at constant p and x a . 

To quantify the degree of overall orientational order¬ 
ing we have calculated the nematic order parameters of 
each species, S^ a \ by diagonalizing the corresponding or¬ 
dering matrix [28] and extracting the largest eigenvalues. 
The corresponding eigenvectors yield the directors n s and 
n r of the DSS and the rods, respectively. From these 
we calculate the degree of magnetization (if present), 


Pi = p'i ' n s \/N s ^ and the biaxiality parame¬ 

ter B = ^ 3 (n r • n s ) 2 — 1 /2^ [29]. Finally, the trans¬ 
lational structure is analyzed by various pair correlation 
functions, see [30] for details. 

We start by considering the case cr* = 1 and x r = 0.8. 
This corresponds to a ratio of volume fractions of </> r /0 s — 
12, which is in the experimentally accessible range [14]. 
As a reference for our GB/DSS simulations we have first 
considered mixtures of GB rods and unpolar soft spheres. 
At the thermodynamic conditions considered, these dis¬ 
play fully miscible isotropic (I), nematic (N) and smectic- 
B (SmB) phases without any indication of global demix¬ 
ing or microphase separation [31] (contrary to, e.g., mix¬ 
tures of GB rods and attractive Lennard-Jones spheres 
[16]). Indeed, the miscibility of our reference system is 
a non-trivial finding given the well-established fact that 
the amount of nanoparticles that can be dispersed in LC 
phases crucially depends on the spheres’ chemical affinity 
[18]- 

We now turn to magnetic mixtures with cr* = 1 and 
x r = 0.8. In Fig. 1(a) we present a corresponding state 
diagram in the T*-p* plane, obtained by MC simulations 
with N = N r + N s = 720 particles (test simulations with 
N = 2000 and N = 4000 particles revealed no apprecia¬ 
ble differences). We note that Fig. 1(a) gives an overview 
of the nature of equilibrated states rather than true phase 
transition lines; obtaining the latter requires extensive 
free energy energy calculations not performed here. 


From the perspective of the LC host matrix, the over¬ 
all state behavior of the system is similar to that of the 
non-magnetic reference mixture. An example of the rod’s 
overall ordering behavior at fixed density is given by the 
plot of the function S ^ (T*) in Fig. 1 (b). On cooling, 
shows a step-like increase from values essentially zero (I 
state) to « 0.9 in the SmB state. We have verified the 
presence of SmB ordering via the pair correlation func¬ 
tion along the rod’s director, g jj r l (ry), which displays 
clear differences between the nematic, translationally dis¬ 
ordered state and the smB state characterized by forma¬ 
tion of layers with local hexagonal-like order (see [30]). 
Quantitatively, we see from Fig. 1(a) that the temper¬ 
atures related to state boundaries are generally larger 
than in the non-magnetic case, with the gap between 



FIG. 1. (Color online) (a) State diagram of a mixture of GB 
rods and dipolar spheres with (7* = 1 (x r = 0.8), involving 
fully miscible isotropic (I), uniaxial nematic (N u ) and uni¬ 
axial smectic-B (SmB) states. The symbols indicate state 
points where the actual simulations were performed. The 
solid and dashed lines indicate state transformations of the 
non-magnetic reference system, (b) Global order parameters 
as functions of T* at p* = 0.40. 



FIG. 2. (Color online) Representative simulation snapshots 
for mixtures with cr^ = 1 and x Y — 0.8 in a) the isotropic 
(I) state [(T*,p*) = (1.1,0.34)], b) uniaxial nematic (N u ) 
state [(T*,p*) = (1.2,0.40)], c) uniaxial smectic-B (SmB) 
state [(T*,p*) = (0.9,0.40)]. In the bottom parts only the 
MNPs are shown for clarity. The directors n r and n s are 
indicated by thick lines. 


the systems increasing with increasing density. In other 
words, the stability of orient at ionally ordered phases is 
enhanced , despite of the absence of an external magnetic 
field. 

To understand the underlying mechanisms we consider 
the microscopic structure of the MNPs in the various LC 
states illustrated in Fig. 2. 


Within the isotropic range [Fig. 2(a)] the dipolar spheres 
tend to self-assemble into worm-like chain structures, yet 
without net alignment [i.e., S ^ « 0, see also Fig. 1(b)]. 
This behavior is familiar from pure, strongly coupled DSS 
systems at low densities [32] (recall that the tempera¬ 
ture range in Fig. 1 corresponds to large dipolar coupling 
strengths of A > 3.5). The strongly correlated nature of 
these chains is also reflected by periodic arcs of the cor¬ 
relation function g^ (r||,rj_) along the dipole moment 
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of a particle, i.e., as function of ry (see [30]). Upon en¬ 
tering the nematic state, the chains spontaneously ” un¬ 
wrap”, yielding essentially straight chains with head-to- 
tail alignment of the dipole moments. Moreover, the axes 
of these magnetic chains form a director n s which coin¬ 
cides (on the average) with the LC director n r . This 
yields an uniaxial nematic (N u ) ordering of the hybrid 
system (B & 1). In fact, the magnetic chains stabilize 
the LC ordering, as indicated by the enhancement of I- 
N u transition temperatures in Fig. 1(a). Vice versa, the 
LC stabilizes the mutual alignment of the chains, which 
is absent in pure DSS fluids at the densities considered 
here [32, 33]. We note, however, that while the chains 
themselves form ferromagnetic entities, the overal mag¬ 
netization is zero, as reflected by the negligible values 
of Pi plotted in Fig. 1(b). Indeed, the absence of mag¬ 
netic correlations beyond a single chain is also reflected 
by the function g\-fi s (rj_), which takes large positive val¬ 
ues near the origin, but then quickly decays to zero (see 

[30] ). Further, there are no indications for a mutual or¬ 
dering of the chains in the plane perpendicular to n r . 
Rather, the chains form a random ” up-down” structure 
[see Fig. 2(b)]. The ordering within the ferromagnetic 
chains becomes even more pronounced upon further cool¬ 
ing the hybrid into the SmB state of the LC host. This 
is directly visible by the snapshots in Fig. 2(c) and also 
reflected by the larger values of the order parameter S^ 
(as compared to the N u state). On the other hand, the 
translational structure of the chains lateral to the direc¬ 
tor remains disordered, contrary to the local hexagonal 
ordering of the rods within the smectic layers [30]. We 
term these novel complex fluids uniaxial ferrosmectics. 

The fact that oriented chains of MNPs inside a ne¬ 
matic LC stabilizes the host’s ordering has been found 
in several experiments [3-6], and has also been theoreti¬ 
cally analyzed by minimizing the elastic free energy [8]. 
The essential idea is that the presence of long magnetic 
chains enhances the strength of the local director field of 
the rods. However, in all these cases the magnetic or¬ 
dering was induced by an external magnetic field. In our 
case, the self-assembly of the MNPs occurs spontaneously 
as a result of the strong dipole-dipole interactions, yield¬ 
ing worm-like chains already in the isotropic regime. As 
soon as these chains have been aligned by the LC host, 
the elastic distortions onto the LC host are similar to 
those in an external field, explaining the LC’s enhanced 
stability. 

Not suprinsingly, the effects reported here are sensitive 
to the size ratio cr* of the spheres relative to the GB rod’s 
width. Indeed, MC simulations (to be reported elsewhere 

[31] ) of systems with significantly smaller dipolar spheres 
(cr* = 0.25) show that the sphere’s structure is character¬ 
ized by rings, chains, and branched structures in all LC 
phases. The nematic ordering of the LC then leads to an 
overall alignment of these highly correlated structures. 
Mixtures with somewhat larger spheres such as cr* = 1.5 



FIG. 3. (Color online) (a) State diagram at size ratio cr* = 2 
(x r — 0.9), involving fully miscible isotropic (I), biaxial ne¬ 
matic (Nb) and biaxial lamellar (Lb) states. The solid line in¬ 
dicates the I-N transformation of the non-magnetic reference 
system. Open symbols indicate defect-rich states (p* = 0.305) 
or hysteresis in cooling-heating series, (b) Functions S^ a \T*) 
at p* — 0.338. Cooling: solid symbols; heating: open sym¬ 
bols. 

behave similar to the case cr* — 1; however, at cr* =2 
the picture changes dramatically. 

As an example we consider systems at x r = 0.9. This 
corresponds to a volume fraction ratio of (j) r /(j) s = 3.4, 
which is at the lower limit achievable in real systems 
[14]. Figures 3(a)-(b) give an overview of the possible 
states (obtained with N = 1251) and the temperature 
dependence of the order parameters, respectively. 


The I state is structurally similar to the case a* = 1 
in that the MNPs form worm-like chains, but no global 
alignment. At low temperatures, both the LC and the 
MNPs exhibit global orientational order, with the char¬ 
acteristics, however, being markedly different from the 
case cr* = 1. First, the boundary between isotropic and 
orientationally ordered states is only slightly enhanced 
relative to the non-magnetic reference case. This already 
indicates that the MNPs exert a stronger perturbation 
on the LC host. The second, more striking difference is 
that the low-temperature phases at cr* = 2 are no longer 
uniaxial, but biaxial in the sense that the director of the 
magnetic chains stands perpendicular to that of the GB 
rods. 

We first consider the biaxial nematic (Nb) state. As 
seen from Figs. 3(b) and 4(a), this state is character¬ 
ized by pronounced nematic ordering of the rods (0.4 < 
S( r ) < 0.7) without any positional ordering along the 
(rod) director. Within this ordered matrix, the DSS ar¬ 
range into chains which tend to align, yet not as pro¬ 
nounced (0.15 < SW < 0.3) as in the case cr* = 1. The 
loose ordering of the chains is also visible from the snap¬ 
shot presented in Fig. 4. The most striking finding, how¬ 
ever, is that the principal directors of the magnetic chains 
and the rods are, on average, perpendicular to each other 
(i.e., B = — 0.5). To check for finite-size effects we have 
performed simulations at N = 2536 with no appreciable 
differences; the directors always arrange spontaneously 
in an orthogonal way. We conclude that the system is 
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in a biaxial nematic state, contrary to mixtures with size 
ratios cr* < 1.5 that exhibit uniaxial nematic states. 

A further new feature not occuring at smaller values 
of <t* concerns the low-temperature behavior. As seen 
from snapshots in Fig. 4, the system forms a biaxial, 
lamellar-like (Lb) state characterized by microphase sep¬ 
aration into alternating, rod-rich and sphere-rich regions, 
with a considerable degree of interdigitation. In detail, 
the rods are arranged in smectic-like layers perpendicu¬ 
lar n r , whereas the dipolar spheres form aligned chains 
in thinner layers arranged in between the smectic layers 
of rods. 



FIG. 4. (Color online) Representative simulation snapshots 
for mixtures with a* = 2 and x r — 0.9 in a) the biaxial 
nematic (Nb) state [(T*,p*) = (1.8,0.338)] and b) the biaxial 
lamellar (Lb) state [(T*,p*) = (1.2,0.338)]. On the right sides 
only the MNPs are shown for clarity. Part b) additionally 
includes a snapshot of the structure in one dipolar layer. The 
directors n r and n s are indicated by thick lines. 


These chains point randomly parallel or antiparallel 
along the sphere director n s _L n r , with the overall de¬ 
gree of alignment being substantially larger (, S^ ~ 0.6) 
than in the nematic phase. The net magnetization is 
zero. Interestingly we have found, by analyzing appro¬ 
priate correlation functions, that the local structure of 
chains within a layer and with neighboring layers some¬ 
what resembles a two-dimensional rhombic lattice in the 
plane perpendicular to the chains’s director [30]. Also, 
we have confirmed the reversibility of the cooling-heating 
procedure by examining the pair correlations [30]; for ex¬ 
ample, both pjj r l (ry) and pjj s l (ry) become structure¬ 
less upon heating from the lamellar state to the nematic 
state. Nevertheless, it should be noted that the trans¬ 
formation is accompanied by considerable hysteresis (see 
solid/open symbols in Fig. 3). 


Lamellar microphase-separated structures have also 
been observed in purely entropic mixtures of hard rods 
and hard spheres [17], as well as in GB-Lennard Jones 
mixtures where the strength of rod-sphere attraction is 
isotropic [16]. In both cases, lamellar phases occur al¬ 
ready at <r* = 1. In the present system, the low- 
temperature state at a* = 1 is not lamellar [see Fig. 2(c)] 
but rather smectic-like with magnetic chains aligned par¬ 
allel to the smectic director. This satisfies not only the in¬ 
teractions within each species but also the rod-sphere in¬ 
teractions which favor planar anchoring. Only at cr* = 2, 
where the spheres perturb the smectic order sufficiently 
strongly, the magnetic particles finally form a layer in 
their own. 

In conclusion, our MC simulation study of binary mix¬ 
tures of GB rods and dipolar soft spheres with compa¬ 
rable size ratios (such as in real colloidal-rod systems 
[13, 14]) suggests intriguing ordering scenarios of these 
model LC-MNP hybrid systems not known so far. In 
particular, we have observed (i) spontaneous formation 
of ferromagnetic chains, yielding an uniaxially ordered 
hybrid without net magnetization but enhanced ordering 
of the LC matrix, and (ii) biaxial nematic and lamellar 
phases, where the magnetic particles arrange into chains 
oriented perpendicular to the LC director and where, at 
sufficiently low temperatures, local translational order¬ 
ing resembling a rhombic structure is formed. This novel 
biaxial LC-MNP system incorporates features of biaxial 
phases in general [34, 35] and, in particular, of biaxial 
mixtures of rod-disk mesogens [29]. Due to the magnetic 
component, both the uniaxial and biaxial ordered hy¬ 
brids can be expected to possess strong anisotropic sen¬ 
sitivity to external magnetic and electric fields, as well 
as to shear forces, with consequences that are yet to be 
explored. Moreover, the different ordering scenarios will 
have influence on light propagation and thermal trans¬ 
port in corresponding real soft magnetic hybrids. Our 
MC results may stimulate further experimental investi¬ 
gations of such systems. 

We thank R. Stannarius and S. Odenbach for stimu¬ 
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This supplementary information is organized in two sections: In Sec. I we present definitions for various pair 
distribution functions which we used to analyze our systems. Numerical results for rod-sphere mixtures with size 
ratios a* = 1 and a* = 2 are presented in Sec. II. 


I. DEFINITIONS OF DISTRIBUTION FUNCTIONS 


To analyze the positional ordering we consider various direction-dependent pair correlation functions in addition to 
the usual radial pair correlation functions (r) [1-3], where a= r (rods) or s (spheres). To start with, we calculate 
the longitudinal correlation function 


*41 (r„) = 


Ei&t(r || ~ lx -n a |) 


( 1 ) 


AV 2 p(N a -l ) 

where AV 2 = tt (r 2 - (rij • ^a) 2 ) Ary is the volume of a cylindical shell with thickness Ary = 0.05cr. The function 

#jj a ^ a (ry) measures the distribution of particles of a specific type along the corresponding director, i.e., ry is the 
projection of the connection vector onto n a . 

Likewise, the distribution perpendicular to the director of one species is measured by the correlation function 


d a k (r±) = ( 


S ( r -L - y r ij - ( r ij ■ *a)' 
AV 3/9 (iV a — 1) 


( 2 ) 


where r± = |— (r^- • n a ) n a |, and AV 3 = 7 tL ^(r_ j_ + A r±) 2 — r 2 ^ with L being the height of the cylindrical shell. 
Further, dipole-dipole correlations along the direction perpendicular to h s are quantified by 


*41 ( r -0 = 


Ei#iM r± ~ A /r 2 . - ( r ir h s ) 


ij U 


(E i#j (n.-) 


(3) 


where cos 6^ — fi i • ftp and fi i is the dipolar unit vector of particle i. 

Finally, to analyze the structure with respect to the dipole moment of one dipolar particle, we calculate the two- 
dimensional correlation function 


g (s) (nvx) = ( 

where AV = tt ((r± + Ar±) 2 — r 2 ^ A ry. 


Ei^i S (x - \fdj - (»*y • Ai) 2 ) S (»-|| - |ry • Ail) 


Al>(.Vs - 1) 


(4) 


II. NUMERICAL RESULTS FOR CORRELATION FUNCTIONS 
A. Mixtures with aj = 1.0 

We start by considering correlations between the rod particles. The structure formation of the rods with decreasing 
temperature at a* = 1.0 (see Figs. 1 and 2 in the main manuscript) is indicated by the correlation function gjj r l (ry) 
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FIG. 1. (Color online) (a) The correlation function g jj r l (ry) and (b) (r_ l) in the uniaxial nematic (N u ) state [(T*,p*) = 

(1.2,0.40)] and the smectic-B (SmB) state [(T*,p*) = (0.9,0.40)] at a £ = 1.0. 
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FIG. 2. (Color online) Two-dimensional pair correlation function (r||,rj_) in the isotropic (I) state at [(T*,p*) = (1.4,0.34)] 
and og = 1.0. 


defined in Eq. (1). In Fig. 1 this function is plotted at density p* = 0.4 for two temperatures pertaining to the nematic 
and the smectic state, respectively. 


In the uniaxial nematic (N u ) state, there are no appreciable correlations. This changes upon decreasing the tem¬ 
perature, where the undamped oscillations in pjj r l (ry) clearly reflect the layering typical of a smectic state. The 

rod correlations (r±) perpendicular to n r [see Fig. 1(b)] indicate smectic-B (SmB) state characterized by local 

hexagonal ordering within the layers. 

We now consider the dipolar particles. In the isotropic (I) state of the LC the dipolar spheres form a (globally 
isotropic) network of wormlike chains (see Fig. 2(a) of the main manuscript). A characteristic signature of this 
arrangement are the pronounced correlations along the dipole vector of a given particle. These are visible in Fig. 2 
where we plot the two-dimensional pair correlation function defined in Eq. (4) at a representative I state. 


The correlations along the chains are reflected by the periodic arcs near the ry-axis. The lack of correlations along 
the r^-axis indicate that neighboring chains are essentially decoupled, as expected due to the small density of DSS. 

Finally, we present in Fig. 3 the influence of temperature on correlations between the dipolar spheres (for snapshots, 
see Fig. 2 in the main manuscript). 
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FIG. 3. (Color online) The dipolar correlation functions (a) g^-fi ( r -0 an d (b) g (rj_) for p* = 0.40 and various temperatures 
T * ranging from the isotropic to the smectic-B regime (cr* = 1.0). 




FIG. 4. (Color online) The distribution functions (a) pjj r l (ry) and (b) gjj s ,l (ry) along the nematic director in the biaxial 
nematic (Nb) state at [(T*,p*) == (1.8,0.338)] and the biaxial lamellar (Lb) state at [(T*,p*) = (1.2,0.338)]. The size ratio is 
<7g = 2.0. The data pertain to the system size N = 2536. 


Figure 3(a) shows that for all temperatures considered, including the lowest ones corresponding to the uniaxial SmB 
state, the perpendicular distribution function (rj_) takes large positive values only near the origin and then 

quickly decays to one (i.e., no correlations). Similarly, the dipole-dipole correlation function (rj_), which is 
sensitive to the particle orientations, is positive (parallel orientation) only at very small distance in perpendicular 
direction. Both functions thus indicate that there are no polar domains beyond a single chain. Moreover, the polar 
chains lack of any long-range translational order. 


B. Mixtures with crj = 2.0 

The mixtures with cr* =2.0 are characterized by the appearance of a biaxial nematic (Nb) and, at even lower 
temperature, a biaxial lamellar (Lb) state (see Figs. 3 and 4 in the main manuscript). To understand the differences 
between these states on the level of correlation functions, we plot in Fig. 4 exemplary results for the rod-rod distribution 
function pjj r l (ry) and the sphere-sphere distribution function pjj s l (ry). 


At T* = 1.8 (Nb state), both distribution functions are essentially structureless, indicating that neither the rods nor 
the spheres are correlated along the nematic director. Upon further cooling into the Lb state, the function (ry) 
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FIG. 5. (Color online) The distribution functions (a) gjj s ^ r (r\\) and (b) g±]ft s ( r ^) th e biaxial lamellar (Lb) state [(T*,p*) = 
(1.2,0.338), o $ = 2.0]. The maxima of the correlation function that are depicted by arrows correspond to distances in the 
sketch of the chain arrangement in Fig. 6. 
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FIG. 6. (Color online) Sketch of the (idealized) arrangement of magnetic chains inside the Lb state (for snapshots, see Fig. 4(b) 
main manuscript). The chains point randomly into or out of the plane. Also indicated is the director of the rods. 


[see Fig. 4(a)] develops oscillations characteristic of the formation of a layered structure. Interestingly, this is also 
reflected by oscillations of the sphere distribution function pjj s l (ry) along the rod n v director [see Fig. 4(b)]; these 
oscillations reflect the formation of dipolar layers between the rods (see Fig. 4(b) in the main manuscript). Indeed, 
the large peak at ry ~ 3.13cr occurs at nearly the same value as that in pjj r l (ry), reflecting that the length of the 
rods is the dominant length scale here. 

Additional information about the distribution of spheres is gained by the function (rj_) which is sensitive 

to structure formation perpendicular to the director characterizing the magnetic chains. Note that this function 
includes contributions from both, neighboring chains within a dipolar layer (see snapshot in Fig. 4(b) (inset) in the 
main manuscript) and neighboring chains in adjacent layers. A plot is shown in Fig. 5(b), while Fig. 5(a) contains 
for comparison the function pjj s l (ry) [same data as in Fig. 4(b)]. 


It is seen that g ( f!f is (rj_) [as does pjj s l (ry)] displays pronounced maxima indicative of the presence of short-ranged 
translational order of the ferromagnetic chains. Analyzing the positions of the maxima of both functions (indicated 
by colored arrows) we can draw a sketch of the arrangement of chains in the Lb state, see Fig. 6. Here, the chains 
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point into (or out of) the plane, whereas the rod’s director lies horizontally. Thus, the chains are positioned within 
layers squeezed between the rods (see snapshot in Fig. 4(b) in the main manuscript). Of course, the sketch illustrates 
a highly idealized situation disregarding thermal fluctuations. 

The spacing between the dipolar layers can be inferred from the positions of maxima in g jj s l (ry) (see violet and 

light blue colors). The second peak of g (rj_) stems from chains in the adjacent layer, on the one hand, and 
from the two neighboring chains inside the layer, on the other hand. Interestingly, this peak (yellow) is positioned 
at a somewhat larger distance (r± « 3.7cr) than that in the longitudinal function (ry « 3.13cr). Thus, the chains 

in adjacent layers are shifted, as indicated in Fig. 6. Finally, the light blue and red peaks in g^f!^ (rj_) stem from 
correlations with chains the next nearest neighboring layer. Altogether, the correlations reflect that the chains are 
locally arranged into a rhombic-like structure. 
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